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Abstract 

In this paper we calculate dipole amplitude and dipole amplitude correlations using 
first and second equation from Balitsky-Kovchegov hierarchy. Our analysis shows that 
even in presence of weak dipole correlation in initial condition mean field approximation 
breaks down through evolution. This difference asymptotically grows not only by absolute 
value but also by relative difference of dipole correlation from mean field value. This could 
affect physical values related to dipole correlation such as elliptic flow V2 and ’back-forward’ 
asymmetry in saturation model. 

1 Introduction 

It is well known that first equation in Balitsky-Kovchegov hierarchy mm in mean field 
approximation have the following form 
_ 2 

dY{T{r,Y)) = ^l dU^2(^^((rN,y)) + (r(r-z,y))-(r(r,K))-(TU,K))(r(r-z,K))) 

( 1 ) 

where {T(r,Y)) is dipole scattering amplitude. 

As was proposed in this equation can be gradually simplified with the following 
transformation 

Tik,Y)= r ^Jo(kr)T{r,Y) ( 2 ) 

Jo r 
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Which lead to simplified form 


d^Y{T{L,Y)) = xi-dL){TiL,Y)) - {T{L,Y)f 


(3) 


Where L = ln{k'^) and integral operator x is dehned by following equation 


x{-dL)T{L) = J^^dL' 


+°° (T{L') - e^-^'T{L) ^ e 


L-L' 


T{L) 


oL-L’\ 


, _ , (4) 

-e— I v'd + e2 (^-^')) 

Solution of equation o was carefully examined by many authors using numerical 
simulation. Moreover as was showed in this equation belongs to Fisher-Kolmogorov- 
Petrovsky-Piscounov (FKPP) equation class [^. Solution of this equation has quite re¬ 
markable properties. For physically acceptable initial condition at asymptotically large Y, 
dipole scattering amplitude {T{k,Y)) can be represented as traveling wave in k space with 
Y as time (hg. ^). Moreover at asymptotically high Y traveling wave form does not depend 
on initial condition. 

However it is not clear if this holds true in case where mean field approximation of initial 
condition is not valid. Since there is no method for analytical solution yet the only way to 
test it is through numerical simulation. Let’s consider first and second Balitsky-Kovchegov 
equations. 


9y(T(r,y)) 


a 



((r(z, y))+ (r(r-z, y)) - (r(r, y)) - (r(z, y)t{t-z, y))) 

(5) 


_ 2 

dY{Tin,Y)T{r2,Y)) = ^l d2^-^^-^^((T(z,y)r(r2,y)) + (r(ri-z,y)r(r2,y))- (6) 

(r(ri,y)r(r2,y)) - (r(z,y)r(ri-z,y)r(r2,y))) 
mruY)T{z, Y)) + (r(ri, y)r(r 2 -z, y)) - 
(r(ri,y)r(r2,y)) - (r(ri,y)r(z,y)r(r2-z,y))), 

Using 0 , simplified form can be obtained 

a^y(T(L,y)) = xi-dL){T{L,Y)) - {T\L,Y)) (7) 

asv(r(Li,y)r(L2,y)) = {x{-dL,)+x{-dLMT{LuY)T{L2,Y)) (8) 

-{t^{Li,y)t{L 2 ,y)) - (r(Li,y)T2(L2,y)) 
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Figure 1: Evolution of dipole scattering amplitude defined by Balitsky-Kovchegov equation. 
log{N{k, Y)) from log{k‘^) dependence for Y = 0,1, 2, 3,4 from left to right 
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This pair of equations allows mean field solution like: 


(r(Li, y) •.. T{Ln, Y)) = A -1 (r(Li, y)) • • • (r(L„, y)) 


Beyond this approximation this equation should be solved numerically. 


(9) 


2 Numerical solution method 


Is is clearly seen that ( 0 © is incomplete since it contain term {T{Li,Y)T(L2, Y)T{L^^ Y)), 
to make this system complete we should make some approximation. Let’s suppose that 


(r(Li,y)r(L2,y)r(L3,y)) = (r(Li,y))(r(L2,y))(r(L3,y)) + (lo) 

((T(Li,y)r(L2,r) - (r(Li,y))(r(L2,y)))(r(L3,y)) 
{{t{L 2 ,y)t{Ls,y) - (r(L 2 ,y))(r(L 3 ,y)))(r(Li,y)) 
((T(L3,r)r(Li,y) - {t{L3,y)){t{l,,y))){t{L2,y)) 

To obtain numerical solution for (00 we should first apply some regularization scheme to 
0 - The most oblivious way for this is cut infinite interval from both sides, ultraviolet and 
infrared. Therefore lets regularize © as: 


xi-dL)nL) = 


dL 


, (T{L') - e^-^'r(L) 


„L-L' 


+ 


T{L) 


( 11 ) 


where << 0 << Lmax- For this regularization it can be shown that residue can 

be safely dropped from numerical computation. 

Lets choose another variable x = , where Ln = ^rnin+Lmax _ Lmin-Lmax 


6L ’ 


Therefore we have for kernel: 


xix)T{x) = J 6 Ldx 


, /T(x') - e^^-^'^^^T{x) e^^-^">^^T{x) 


I 2 _ q{x—x')SL I 


+ 


_|_ q2{x-x')5L 


( 12 ) 


Defining yi(x,y) = (T(x, T)) and iV2(xi, X2, T) = (r(a:i, y)r(x2, y))-iVi(a:i, y)iV2(a:2, T) 
and using m equation (0 can be rewritten as 


d^YNi{x,Y) = x{x)m{x,Y) - NUx,Y) - N 2 {x,x,Y) ( 13 ) 

daYN2{xi,X2,Y)) = (xixi) + xix2))N2{xi,X2,Y) ( 14 ) 

-2N2 (xi, X2, y) (xi, y) - 2N2 {xi, X2, Y)m ix2, Y) 
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Now there is oblivious way to solve this equation. Lets approximate Ni{x,Y) and 
N 2 {xi,X 2 ,Y) with 


Nr, 


Niix,Y)= Y. CliYnix) 


0 

Nmax Nmax 

N2{xuX2,Y)= E E CYY)Uiix^)Ujix2) 

i=0 j=0 

where Ui{x) is some set of orthogonal functions with 

J ^ Ui{x)Uj{x)^{x)dx = 6ijHi 
With this equation (ITTOl) can be rewritten as 


(15) 

(16) 

(17) 


1 


1 


davCHY) = -XijCiiY) - -FijkCtiY)) 


(18) 


d^yC^Y) = + -XjkCfiY) - (19) 

Hi Hj 

2 2 

yii yikl T? /^kl jp 

- ^klj - ^kli 

Hi Hj 

where 

Xij = J ^ Ui{x)x{x)U j{x)H{x)dx (20) 

Fijk = J ^ Ui(x)Uj{x)Uk(x)H{x)dx (21) 

System of ordinary differential equation mm can easily be solved by Runge-Kutta- 
Fehlberg method. 

The most oblivious choose for orthogonal basis Ui(x) is the first kind Chebyshev polyno¬ 
mials. However in this case solution is unstable. The instability comes from points |x| = 1. 
Same applies for second kind Chebyshev polynomials. Therefore we choose Legendre poly¬ 
nomials as orthogonal basis for our purposes. 


3 Results 

For numerical computation we set Nmax = 64, Lmm = —10, L^ax = 50, = 0.2. 

p-k/Qs-ik/Qa)^ 1 

Niik,0) = -(22) 

N2(ki,k2,0) = eNi{ki,Y)Ni(k2,Y)exp{-log'^(ki/k2)) 
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Figure 2: Evolution of dipole correlation defined by Balitsky-Kovchegov equation. N 2 {ki, ^ 2 , Y) 
from log{kl), log^k^) dependence for Y = 0,1, 2, 3,4, 5 from up to down 

with e = 10“^. Numerical solution for absolute and relative correlation N 2 {ki, k 2 ,Y) is 
shown on fig. 2 and fig. 3 respectively. 

It is clearly seen that {T{ki,Y)T{k 2 ,Y)) is smaller what value {T{ki,Y)){T{k 2 ,Y)) 
supposed by mean field approximation. Moreover not only absolute difference between 
{T{ki,Y)T{k 2 ,Y)) and {T{ki,Y)){T{k 2 ,Y)) does not vanish with Y but so does relative 
difference. This difference affect physical values related to dipole correlation such as elliptic 
flow V 2 and ’back-forward’ asymmetry in saturation model. Recently Balitsky-Kovchegov 
hierarchy was reformulated [2] using stochastic field theory. It was showed what m should 
contain additional ’border’ term which related to fluctuation. It therefore possible that 
N 2 {ki, k 2 ,Y) behaviour will be changed drastically. This required further investigation. 
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Figure 3: Evolution of dipole correlation defined by Balitsky-Kovchegov equation. 

^ 2 (^ 1 , ^ 2 , y)/Ni{ki^ Y)Ni{k 2 -, E) +1 from log{k‘f), log{kl) dependence for Y = 0,1, 2, 3,4, 5 from 


up to down 
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